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D COMPLEX*16 Sine and Cosine 


HISTORY ; Detailea Current Edit History 
DECLARATIONS 

MTHSCDSIN = D COMPLEX*16 SINE 

MTHSCOCOS = D COMPLEX*16 Cosine 

WORKER = do all the work 
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MTHSCOSINCOS D COMPLEX*16 Sine and Cosine 16-SEP-1984 01:07: AX/VMS Macro V04-00 Page 
1-005 at 94:90 38 YRTHRTL. SREIRTHCDSINC MAR: 1 ° 
TITLE ye aaa D COMPLEX*16 Sine and Cosine 

IDENT /1-002/ 3 File: MTHCDSINC.MAR Edit: RNH1002 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
0 MS OF CH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
SORPORATION NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR _ RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


eaeeeeeneeeeeeeeaneneaenannne 


PRRRARASASAAAAAALALALAALASASARSAEALELASASA ARSE ALAS ASESES RRR ARRAS SESSA RAD SS OY 


FACILITY: MATH LIBRARY 


ABSTRACT: 
Return the SINE of 


A OMPLEX*16 number 
Return the COSINE of 


dc 
a D COMPLEX*16 number 


VERSION: 1 
HISTORY: 


AUTHOR: 
Steven B. Lionel, 26-July-1979 


MODIFIED BY: 
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-SBTTL HISTORY ; Detailed Current Edit History 


tis 


09 1-00) = Adapted from HTHSCSINCOS version 1-002, SBI 26-July=1979 


l= 


2 = Changed shared external references to G6* R 


BTWSCDSINCOS 6 Sine and Cosine ms 1 ~$E0-1984 91 :07 ‘23 YAX/WRS Ma Macro V04-00 Page by es 


d¢ *] 
DECLA ONS ~SEP=1984 MTHRTL.SRC SMTHCDSINC.MAR; 1 
-SBTTL DECLARATIONS 


; INCLUDE FILES: 


suisse SYMBOLS: 
ABL GBL 


MACROS: 


; PSECT DECLA pn S: 
~-PSECT _MTHSCODE PIC, SHR, LONG, EXE, NOWRT 


EQUATED SYMBOLS: 


C009 09 090909 NINN NINN NINO AAA AOAAOOOwn 


MEW SO OONAUE WN" OOONAOULS WN —O”o 
eee 


> OWN STORAGE: 
NONE 
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MTHSCOSINCOS D COMPLEX*16 Sine and Cosine 16-SEP-1984 01:07: AX/VMS Macro V04-00 Page 4 
Sr O0S MTHS CDS iN - D COMPLEX*16 SINE 66-SEP-1984 01:90 ike EATHRTL . SRESATHCDSINC MARS 1 . (4) 
6 i -SBTTL MTHSCDSIN = D COMPLEX*16 SINE 
p++ 
0 ys : FUNCTIONAL DESCRIPTION: 
o 9¢ : MTHSCDSIN computes the SINE of a D COMPLEX*16 number (r, i) as 
7) 6 : 
000 Be : result = (SIN(r) * COSH(i), COS(r) * SINH(i)) 
44! 96 ; CALLING SEQUENCE: 
38 $3 3 CALL MTHSCOSIN (result.wdc.r, arg.rdc.r) 
99 ; 
000 100 ; INPUT PARAMETERS: 
00000008 4 +4 arg = 8 ; D COMPLEX*16 argument by reference 
44 108 3 IMPLICIT INPUTS: 
000 104; NONE 
0000 105; 
0000 106 ; OUTPUT PARAMETERS: 
464 107 ; 
00000004 +443 138 result = 4 ; D COMPLEX*16 result by reference 
0000 110 ; IMPLICIT OUTPUTS: 
0000 111; NONE 
0000 a 3 
0000 1135 ; COMPLETION CODES: 
0000 114; NONE 
0000 115; 
0000 116 ; SIDE EFFECTS: 
0000 117; Signals: Reserved Operand if r or i are invalid (-0.0) 
0000 118; MTH$_SINSIGLOS if iri > 2ePI#2e*31, 
0000 119; Floating Overflow if i > 88.028. 
44 120 ; 
000 121 ;-- 
0000 \$¢ 
0000 i 
OOFC 0000 124 -ENTRY MTHSCDSIN, “M<R2,R3,R4,R5,R6,R7> 
00000039'EF 16 0008 125 JSB WORKER : RO-RI ="SIN(r) 
000 1 § 3 R2-R3 = COS(r) 
0008 1 > R4-R5 = SINH(i) 
44 128 3: R6-R7 = COSH(i) 
50 56 64 000 129 MULD R6, RO 3 a) = SIN(r) * COSH(i) 
52. 54 64 43 130 MULD R4, R2 3 R2-R3 = COS(r) * SINH(i) 
54 04 AC 0 00 131 OVL result(AP), R4 3 get result address 
84 50 D 001 1 ¢ mMOVQ RO, (R4)+ ; Store real part 
64 52 7D 001 1 MOVa R2, (R4) ; Store imaginary part 
04 0018 134 ET 
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COMPLEX*16 Cosine 6-SEP=19 MTHRTL.SRCJMTHCDSINC.MAR; 1 
$ -SBTTL MTHSCDCOS - D COMPLEX*16 Cosine 


ov 


p++ 

§ 3 FUNCTIONAL DESCRIPTION: 

MTHSCDCOS computes the COSINE of D COMPLEX*16 number (r, i) as follows: 
result = (COS(r) * COSH(i), -SIN(r) * SINH(-i)) 


CALLING SEQUENCE: 
CALL MTHSCOCOS (result.wdc.r, arg.rdc.r) 


INPUT PARAMETERS: 
arg = 8 


00000008 ; ; D COMPLEX*16 argument by reference 
> IMPLICIT INPUTS: 
3 NONE 
: OUTPUT PARAMETERS: 
00000004 : result = 4 ; D COMPLEX*16 result by reference 
IMPLICIT OUTPUTS: 
NONE 


COMPLETION CODES: 
NONE 


SIDE EFFECTS: 
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CODD OI BODOG SII NN DD DS DD DDD DS TT BB SE EEA 
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Signals: Reserved Operand if r or i are invalid (-0.0) 
MTH$_SINSIGLOS if ir! > 2*PI*2e*31, 
Floating Overflow if i > 88.028. 
OOF C -ENTRY MTHSCDCOS, “M<R2,R3,R4,R5,R6,R7> 
00000039°EF 16 JSB WORKER : RO-R1 ="SIN(r) 
1 : R2-R3 = COS(r) 
1 3 R4-R5 = SINH(i) 
1 ; R6-R7 = COSH(i) 
56 52 64 1 MULD2 R2, R6 3 R6-R7 = COS(r) * COSH(i) 
2? 28 re 4 MNEG RO, RO : RO-R1 = -SIN(r) ; 
52 4 § 7 MULD RO, R4, R2 ¢ R2-R3 = -SIN(r) * SINH(i) 
50 56 4 bp 8 mMOVa R6, 3: RO-R1 = COS(r) * COSH(i) 
54 04 AC 0 MOVL result(AP), R4 ; get result address 
84 23 D 80 MOVQ RO, (R4)+ ; Store real part 
64 2 7d 003 MOVQ R2, (R4) ; Store imaginary part 
04 0038 RET 
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1-00 


6D 00000000'GF 


50 08 AC 
50. = 08 AD 
00000000 * GF 
52 08 50 
54 50 52 
7E 54 00 
54 50 52 
7E 54 00 
50. ~— 08 BC 
00000000 ' GF 
7E 50 

50 =—08 BC 
00000000 ' GF 
52s BE 

56 Ss BE 
548 
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-SBTTL WORKER = do all the work 


: Setup eer handler 
; Com mpute HF 


zw 
= 
‘6 
z2wD 
NVI 
nun 


WORKER: 


MTHSFLAG_JACKET 
MOVAB G*‘MTHSSJACKET_HND, (FP) 


arg (AP), RO 
MOVD BCRO) 

JSB GoM THSDEXP. R7 
DIVD3 RO, #1, R2~ 


SUBD3 Ree a, — 

#0.5, R4, =(SP) 
ADDDS =R2, RO, R4 
MULD3 0.5, R4, -(SP) 
MOVD @arg(AP), RO 
JSB GoM HsDcOs R7 
MOVQ RO, -(SP) 


MOVa Borate? 
J Husosini 87 


MOVQ (SP)+, R6 
MOVQ (SP)+, R4 
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; set up error handler 


; set handler address to jacket 


handler 


+ EXP(-i) 
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MTHSCOSINCOS D COMPLEX*16 Sine and Cosine 16-SEP-1984 01:07:23 VAX/VMS M v04-00 P 7 
Symbol table grSepiobe 91220588 EMTARTe eeeiarucoseuc.maRs1 2% ¢S 
ARG = 00000008 
MTHSSJACKET_HND teeeeeee =X (1 
MTHSCOCOS $08 019 RG a 
MTHSCOSIN 0000000 RG 1 
MTHSOCOS_R7 eterereer =X 00 
MTHSDEXP_R7 tererene x 09 
MTHSDSIN_R7 eereeeee 8 8X 0 
RESULT = at 
WORKER 00000039 R 01 

teaser sens en enn + 

: Psect synopsis ! 

pe oc meweseecceee= + 
PSECT name Allocation PSECT No. Attributes 
s 0S ; 0000000 ( 0.) 00 ¢ OO.) NOPIC USR CON’ ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 0000008 c teeek Oh 4 Bod PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 

Peers nme m neem eee acee naan + 
; Performance indicators ! 

Phase Page faults CPU Time Elapsed Time 
Initialization 33 00:00:00.09 00:00:01 -66 
Command processing 111 eS Ba 00:00:03.38 
Pass 1 89 00:00:00.7 00:00:02.56 
Symbol table sort 0 Se $4 00:00:00.01 
Pass 2 54 0:00:00.61 00:00:03.07 
Symbol table output § 0:00:00.03 00:00:00.25 
Psect synopsis output 0:00:00.02 00:00:00.08 
Cross-reference output ° 00: 3 :39-9? Bo 88 08 «28 
Assembler run totals 29 00:00:02.14 0:00:10.62 


The working set Limit was 900 pages. 

3610 bytes (8 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 9 non-local and 0 local symbols. 
source Lines were read in Pass 1, producing 14 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 


See merece er eee new memo we} 


Macro library name Macros defined 
“S25SSDUAZB:(SYSLIBISTARLET.MLB;20 0 SS 0 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHCDSINC/OBJ=OBJ$:MTHCDSINC MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MS 
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